Quantum phase transition of Bose-Einstein condensates on a ring nonlinear lattice 
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We study the phase transitions in a one dimensional Bose-Einstein condensate on a ring whose 
atomic scattering length is modulated periodically along the ring. By using a modified Bogoliubov 
method to treat such a nonlinear lattice in the mean held approximation, we find that the phase 
transitions are of different orders when the modulation period is 2 and greater than 2. We further 
perform a full quantum mechanical treatment based on the time-evolving block decimation algorithm 
which confirms the mean field results and reveals interesting quantum behavior of the system. Our 
studies yield important knowledge of competing mechanisms behind the phase transitions and the 
quantum nature of this system . 
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I. INTRODUCTION 



In the past few years, ultracold atoms confined in optical lattices have generated a great amount of excitement 
in the physics community. They provide the unique opportunity to realize various many-body models that are of 
fundamental importance in physics More recently, nonlinear lattices formed by periodically modulating atomic 
interaction strengths have also received much attention. A comprehensive review of nonlinear wave phenomena 
supported by nonlinear lattices can be found in Ref. dj. There are two basic physical systems that can potentially 
I ' realize nonlinear lattices, both of which can be described by nonlinear Schodinger equations. One uses electromagnetic 
waves subject to inhomogeneous nonlinear optical media. The other one is based on atomic Bose-Einstein condensates 
(BECs) with modulated s-wave scattering length. In this work, we will focus on the latter system although much of 
the physics are common to both. 

We extend ourprevious work reported in Ref. [1| to study a BEC on a one-dimensional (ID) nonlinear ring-shaped 
lattice. In Ref. [3[, we considered a BEC on this ring lattice whose atomic scattering length is modulated according 
to 



a(0) = a o sin(d0), (1) 



where is the azimuthal angle along the ring and d — 2 is the spatial modulation frequency. We have shown that, as 
• , the modulation depth clq is increased, the condensate can undergo a second-order symmtry-breaking quantum phase 
transition from a soliton-like state to a spatially periodic condensate that matches the scattering length modulation. 
In the present work, we generalize our investigation to a larger spatial modulation frequency d > 3 and compare the 
results to the d = 2 case. We developed a new mean-field technique to study the semi-classical behavior of the system, 
• • ' and found that a similar symmetry-breaking phase transition occurs for d > 3 as the modulation depth is increased. 
However, the phase transition is now of first order. We also carried out a numerical full quantum mechanical treatment 
of the system based on the time-evolving block decimation (TEBD) algorithm. Both static and dynamical properties 
, of the system are investigated. 

a ■ 

II. THE MODEL HAMILTONIAN 

The system considered here is similar to that studied in Ref. 0, 0] ■ N bosons are confined in a toroid of radius R 
and cross sectional area S. By sufficiently tightening the radial confinement and freezing the atoms in that direction, 



'Electronic address: zwzhou@ustc.edu.cn 
^Electronic address: xizhou@ustc.edu.cn 
t Electronic address: hpu@rice.edu 



2 



we can treat the atoms as a one dimension system on a ring. The Hamiltonian of the system can be written in the 
following dimensionless form: 



H = 



d0 



u 



^ (6>) ^ (0) ip (0) ip (0) 



(2) 



where the first term in the integral represents the kinetic energy and the second the interaction energy. For simplicity 
in notations, we measure energy in units of H 2 /(2mR 2 ). The dimensionless interaction energy U (0) = 8ira (0) R/S, 
where a (0) is the periodically modulated s-wave scattering length. In our work, we consider the situation where the 
scattering length between atoms is modulated along the ring with d periods: a (0) = ao sin (dO). As in Ref. 0], we 
define the dimensionless interaction strength as: 



7 (0) = -^^ = 7osin (rf6-), 



(3) 



where 70 = —Aa^RN / S represents the modulation depth of the interaction parameter 7 (0). 

By taking the Fourier expansion of the field operator ip (0) = J2k ^^ eik6 ® k > wnere k takes integer values in order 
to satisfy the periodic boundary condition and at is the bosonic annihilation operator for plane wave mode with 
wavenumber fc, the Hamiltonian can be rewritten as: 



H = k 2 a\a k + ^ a\a\a m an{5d+ m +n-k- 

k klmn 
—$—d+m+n-k-l)- 



Since the number of atoms is fixed, we have N = ajai. Because of this, the kinetic energy term in the Hamiltonian 
can also be written as jfj^ki k 2 a\aka\ai. 



III. MEAN-FIELD TREATMENT 



In this section, we first consider the mean- field solution valid for iV ^> 1. In this case, the kinetic energy term can 
be approximated as 

— k 2 a\a k a\ai « — ^ k 2 a\a\a k ai , 

k,l k,l 

and the total Hamiltonian can thus be cast into a biquadratic form as: 

H = — ^ a ijkl alapkai . (4) 

i,j,k,l 

We will now describe a modified Bogoliubov method we use to find the stationary solution and the excitations of the 
system. 



A. Modified Bogoliubov Approach 

In the absence of atomic interactions, the ground state is quite trivial: all the atoms occupy the zero-momentum 
mode eto. In the presence of interaction, this is no longer true. However, we may conjecture that the system condenses 
into a different ground mode \o- This mode \o, together with other orthogonal modes {xi}'s that form a complete 
set, are related to the {a^} modes through a unitary transformation U: 

(Xo, Xi,X2, ■~) T = U (a , ai, a 2 , ...) T . (5) 
In terms of {Xi} & n d {xlli the biquadratic Hamiltonian in Eq. takes the following form: 
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H = X0X0X0X0 + 2^ c kiXkXi + h.c. + ^ d k iX k Xl 

y fc,;#o y fcj^o 

+ ( ^ X! PfcimxIxiXm + /l.C. J + qklmnxlxlXmXn + 

y k,l,m^O y k,l,m,n^Q k,l 



QklmnXkXlXrnXn + ^ 2^ nwXfcXi ■ (6) 
fe,i 

This Hamiltonian can be simplified by a few considerations. First, it is assumed that most atoms will be in the 
condensate mode xo- Under the mean field approximation, operators for the macroscopically occupied condensate 
mode are replaced by c- numbers, i.e., Xo> Xo ~^ \fN ■ Since occupation numbers in other \k modes are very small, we 
can drop terms involving 3 or more operators in Xk and x\ f° r ^ 7^ 0. After this exercise, we obtain the following 
effective Hamiltonian up to second order in {xk} and {xll : 

H eS = c N + ^2 c kiXkXi + h.c. + ^2 dkiXkXl ■ ( 7 ) 

\k,l^0 J k,ljtO 

This effective Hamiltonian can be diagonalized by the Bogoliubov transformation and the system's elementary 
excitaitons are quasiparticles in nature. In order to investigate the stability of the system, we need to analyze 
the energy spectrum of these quasiparticle excitations. For this purpose, we should work with the following grand 
canonical operator to acccount for the conservation of atom numbers: 

K = H eS -fiN. (8) 

Here, [i is the chemical potential and can be calculated from the condensate energy E = (H) w cqN as: 

8E dc , T 

^9N = C0 + M N - (9) 

Now, we may diagonalize the operator K by using the Bogoliubov transformation on {xfc} and {xl} an d obtain 
the excitation spectrum of quasiparticles. (This will be elaborated on later.) If there is no imaginary excitation 
frequencies, we claim that the condensate mode xo is dynamically stable. 

In the above prescription, the key step is to search for the appropriate unitary transformation defined in Eq. (J5j) that 
transforms Hamiltonian Q into ©. Although there are a great number of unknown parameters in the undetermined 
unitary matrix U, further analysis shows that only the elements in the first row of U are necessary for determining 
the form of the Hamiltonian ©. 

To see this, we note that in order to transform Hamiltonian Q into ^ via the unitary matrix U, a fundamental 

requirement is to maintain the biquadratic terms of the operators ^xj, Xo^ and to eliminate all the cubic terms. To 
this end, we may use the operators (x\>Xj^ to represent the operators faj,a m J in Hamiltonian ([?]) via: 

a i = ^ ] u ji Xj 

j 

where tt,-j is the matrix element of the unitary matrix U , and we obtain the following two equations : 

CO = ai l kl U 0i U 0j U 0k U 0l > ( 10 ) 

°= Y (a ljk i + aijik) [u 0l u 0j UQ k (ai - u^xo)} ■ (11) 

i,j,k,l 

Since xo = J2i u oi a i, Eq. (fTT|) can be recast into a set of equations in operators {ai} (the number of this set of 
equations depends on the cut-off of the Bose modes) and can be solved by numerical method. Once the representation 
of the operator xo is determined, the parameter Co and the chemical potential /i can be obtained by solving Eqs. (|10[) 
and ©, respectively. Therefore, Eqs. (|10l) and (|TT|) together represent an algebraic form of the Gross-Pitaevskii (GP) 
equation. In Appendix [XI we will show that they are indeed equivalent to the ordinary GP equation. The main 
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advantage of Eqs. (fT0|) and (fTTj) is that, in principle, all the stationary states (both dynamically stable and unstable 
ones) of the system can be found. When more than one stable solutions are found, the one that is dynamically 
stable and with the lowest energy will be identified as the ground state of the system. In contrast, with ordinary 
GP equation, using imaginary time evolution method one can only find the dynamically stable states, and often just 
the ground state. Therefore, Eqs. (jlOp and (JTTJ) are supeior for studying the phase transitions in our system, where 
information beyond the ground state is needed. 

Once the condensate state is determined, the Bogoliubov spectrum of quasiparticle excitations above it can be found 
by diagonalizing K defined in Eq. © using the Bogoliubov transformation. The details are described in Appendix iBl 

B. Mean-field quantum phase transition 

Using the modified Bogoliubov method outlined in the previous section, in principle we can find all the stationary 
states (not just the ground state) and their excitation spectrum for any modulated atomic interactions. This provides 
a more thorough picture of the energy landscape of the system and deeper insights into possbile quantum phase 
transitions induced when certain parameters are varied (in our case, the modulation depth of the interaction strength 
7o)- 

Our goal is to study the mean field quantum phase transition for different modulation period d as the modulation 
depth 70 is varied. We concentrate on the low-energy states, meaning stationary states with energy close to that of 
the ground state. We find that there are mainly two types of stationary states in the low energy regime. One type 
has a density profile matching the modulated period of the scattering length. We refer to such states as symmetric 
states. The other type features a density profile that spontaneously breaks the symmetry of the modulation. We 
refer to such states as asymmetric states. The asymmetric states is always (i-fold degenerate with the peak density 
located at ^-^j {i = 0, 1, ...d — 1), where the local interaction energy U{9) reaches the minimum. Under the mean 
field treatment (MFT) , the symmetry-breaking quantum phase transition from the symmetric type to the asymmetric 
type have been studied for uniform scattering length [|| and for d — 2 periodic scattering length Q. Here, by taking 
advantage of the modified Bogoliubov method, we investigate the critical point of the quantum phase transitions and 
the prime mechanism driving such quantum phase transitions for arbitrary modulation period d. We find that there 
is a fundamental difference between d = 2 and d > 2. 



o 




0.3 0.4 0.5 0.6 0.7 0.8 

, , , , , , , , , , , 




0.3 0.4 0.5 0.6 0.7 0.8 

To 



FIG. 1: (Color online) Upper panel: Energy of the condensate for d = 2. The red solid line is for the dynamically stable 
symmetric state, the red dashed line is for the dynamically unstable symmetric state and the blue solid line is for the stable 
asymmetric state. Lower panel: Chemical potential of the ground state. At the critical point 70 = 0.528 (indicated by arrows in 
the plots) the ground state changes from the symmetric to the asymmetric type. Shown in the insets are typical wave functions 
for symmetric and asymmetric states. 
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Figure Q] shows the energies and chemcial potentials of low-energy Bose condensate states by varying the parameter 
7o for d — 2. For small 70, the ground state is a symmetric state. As 70 is increased, a symmetry-breaking phase 
transition occurs at a critical value of 70 = 0.528. At this point, the symmetric state becomes dynamically unstable 
and the ground state changes to an asymmetric state. The ground state chemcial potential shows a kink at this 
critical point, wheras the ground state energy curve is smooth. This represents a second-order phase transition Q. 
A similar behavior is also found in attractive BEC with unmodulated scattering length 
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FIG. 2: (Color online) Same plots as in Fig. [T] for d — 3. The symmetry breaking phase transition occurs at the critical point 
70 = 0.85 (indicated by the arrows in the plots) and the dynamical instability sets in for the symmetric state at 70 = 1.04. 

Next, we turn our attention to the case of d — 3. The energies and chemical potentials as functions of 70 are plotted 
in Fig. [5] Similar to the d = 2 case, for small 70, the ground state is symmetric. As 70 is increased, a symmetry- 
breaking phase transition occurs at a critical value of 70 = 0.85. However, unlike in the d — 2 case, at this critical 
point, the symmetric state remains dynamically stable although the energy of the asymmetric state drops below that 
of the symmetric state. Furthermore, the ground state energy curve as a function of 70 shows a kink at this critical 
value, while the ground state chemical potential becomes discontinuous at this point. Therefore, the phase transition 
at this point is of first order. The dynamical instability of the symmetric state does not occur until 70 = 1.04. 

To summarize our mean-field results, we have found that the nature of the phase transitions changes for different 
modulation period d due to the presence of competing mechanisms in this system. When the modulation period d = 2, 
the symmetry-breaking phase transition is of second order and is induced by dynamical instability of the associated 
states. For d = 3, in contrast, the symmetric-breaking phase transition is of first order and is driven by the level 
crossing of different types of states. We have also investigated the cases for d = 4 and 5 and found similar behavior 
as in d = 3. 





IV. QUANTUM MECHANICAL TREATMENT 

So far, we have limited our discussion to the mean-field approximation. Now we perform a fully quantum mechanical 
examination of the system using a numerical method. In our previous work [3j, exact diagonalization is used for 
this purpose. Here, we will use the TEBD algorithm [Bj]. This has a two- fold advantage compared to the exact 
diagonalization method: (1) It allows us to treat larger systems; and (2) in addition to the static properties, we can 
also use the TEBD method to study the dynamical behavior of the system. 

We first discretize the space by introducing an equidistant grid 9i = iA9, (i = 0,1, ...M — 1). We then replace 
the field operator ip (6i) by ip%/y/ A6>, where ipi is a bosonic annihilation operator. In doing so, integrals can be 
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replaced by sums and the second derivative in the kinetic energy term can be approximated by the difference quotient 



ip (Si+i) + "0 — 2^ (9i) /A8 2 . Finally, the discretized Hamiltonian is 



Af-l 



M— 1 



M— 1 



* = E + + ^ E « + E ^ 



A6> 2 



2A6> 



(12) 



Here, the periodic boundary condition leads to the relation ipo = V'm- For TEBD algorithm with periodic boundary 
condition, we refer to Ref. [7|. In our numerical treatment, we typically divide the ring into M = 60 equidistant grids. 
Our code is adapted from the open source package maintained by the group of Lincoln Carr [8] . 



A. many-body ground-state energy 



In Fig. [3J we plot the ground-state energy per atom of the many-body system as functions of 70 for the modulation 
period d — 2 and d — 3 . We see that the ground-state energy curve approaches the MFT result as N increases, which 
is consistent with the usual quantum-semiclassical crossover behavior for finite-size quantum systems. 
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FIG. 3: (Color online) Many-body ground-state energy per atom for iV = 6, 12 and modulation period d = 2 (upper panel) and 
d — 3 (lower panel). The black solid lines are the mean-field results. 



d = 3 




B. quantum correlation 



In the quantum mechanical treatment, unlike the mean field results, the spontaneous symmetry breaking of density 
distribution of ground state wavefunction in real space dose not occur. The density profile of the quantum mechanical 
ground state always matches the spatial modulation of the scattering length Q. However, we can still gain important 
insights into the change in the characteristics of the wavefunctions by examining the quantum correlation which is 
neglected in the mean-field study. 

For further discussion, we define the partial number operator between the interval 9 € [ifi, tpf] as: 

1 l" fl — ~ 

*W/] = 2^y MtfW)1>{6). (13) 
For the particular case of d = 3, we define three partial particle number operators as follows: 

«< = «[(<_!) "p,^], *=1,2,3. (14) 



7 




Using these number operators, we can define the bipartite and tripartite correlation functions as: 

(njfij) 
(n,-) (rij) ' 

(nin 2 n 3 ) 
(m) (n 2 ) (n 3 ) ' 

We plot the bipartite and tripartite correlations as functions of 70 in Fig. [4j All these correlations are monotonically 

(2) (2) (2) 

decreasing functions of 70. In our calculation, the three two-body correlation functions g 12 , S13 and are essentially 
identical, which is also expected from the symmetry of the system. At 70 = 0, the ground state is exactly known: 

( a n N 

I*) ground = l vac )- The theoretical values of bipartite and tripartite correlations can be obtained as: 

« 00 1 
4?(7o=0) = l-^£^ 

1=1 

p <3, (7.=o) = i-^f;i + o(i 5 ), 

1=1 v 7 

which are in good agreement with the numerical results. When N goes to infinity, all the bipartite and tripartite 
correlations approach unity at 70 = 0. 

Figure 0] shows that although both bipartite and tripartite correlations decay as the interaction parameter 70 
increases, the tirpartitle correlation g^ decays into zero much faster than the bipartitle correlation gf^ ■ For the case 

of N = 6 as illustrated in Fig. 01 g^ is essentially zero at 70 = 1.6, while all the gjy are clearly non-zero at the same 
7o- This is reminiscent of the three-body entangled VF-state, which can be written as 

\W) = i= (|100) + |010) + |001}) . 

For the VF-state, the tripartite entanglement characterized by the 3-tangle disappears and the bipartite entanglement 
characterized by concurrence remains finite [9|. For large 70, our mean-field treatment presented earlier reveals that 
the ground state is characterized by asysmmetric state with three- fold degeneracy. Each of these degenerate mean- field 
state features a density peak at 9 = (1 + iin)/6 (i = 1, 2, 3). In the quantum treatment, this degeneracy is lifted by 
quantum fluctuations, and the non-degenerate quantum ground state may be regarded as roughly a W-state formed 
by these three mean-field states. 



(2) 

9ij 



7 (3) 



C. Single-particle density matrix 



Another important quantity to characterize the many-body state is the single-particle density matrix p^ 1 ' whose 
matrix element is defined as [HI [ll| : 



= (®&) . (15) 
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where the expectation value is calculated with respect to the ground state obtained using the TEBD method. Roughly 
speaking, represents the probability amplitude of finding one particle at site i and at the same time another particle 
at site j. 
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FIG. 5: (Color online) The largest three eigenvalues of the single-particle density matrix p' 1 ' for N = 6. 



Because the matrix p^ is Hermitian it can be diagonalized as: 

pW=Nj2PiV*<Pi, (16) 

i 

where the eigenvalues pi are non- negative and satisfy the constraint YliPi ~ 1- F° r a bosonic system with TV ^> 1, 
Pi is closely related to the condensate fraction of the system. It is easy to see that, the system will be in a simple 
condensate state if and only if the largest eigenvaule is of order unity and all the other p^s are of order 0(N~ 1 ). If 
there are multiple p^s of order unity, the system is said to be in a fragmented condensate because all the corresponding 
states have appreciable occupation numbers. Finally, if all the p^s are of order 0(N~ 1 ), then the system is not Bose 
condensed. It is reasonable to speculate that, at small TV, the condensate fraction versus the strength parameter 70 
exhibits quantum crossover behavior. For N — 6 with the modulation periods d = 2 and d = 3, we plot the three 
largest p^s versus 70 in Fig. [5j The figure shows that at small 70, the system may be characterized as a simple 
condensate. It becomes more and more fragmented as 70 is increased. In the large 70 limit, the eigenvalues approach 
some steady state values and there are d eigenvalues which are much larger than the rest. This is consistent with the 
earlier argument that at large 70, the quantum many-body ground state can be roughly regarded as superpositions 
of the ci-fold degenerate mean-field ground states. 

To further quantify the crossover from a simple condensate to a fragmented condensate, we adopt two methods 
as described below. In Fig. HJa), we plot dpo/^7o as a function of 70, where po is the largest single-particle density 
matrix eigenvalue. In Fig. HJb), we plot the overlap of the ground-state wave function pjj |(G(7o)|G(7o + <$7))| for 
a small value of £7 = 0.02, where |G(7o)) denotes the ground state at 70. Both of these quantities measure how fast 
the characteristics of the ground state change as 70 is varied. These two measures provide consistent results: both 
exhibit a dip at some critical value 70 = 0.84 for d = 2 and 70 = 1.08 for d = 3, which we may define as the critical 
modulation depth where the system crosses from a simple condensate to a fragmented condensate. 
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FIG. 6: (Color online) (a) Derivative of the largest single-particle density matrix eigenvalue with respect to the interaction 
strength modulation depth and (b) overlap of the ground-state wave function versus 70 for the particle number N = 6. 



D. Time evolution of the survival probability 



Up to now, we have focused our attention on the static properties of the ground state. The low-energy excitations 
of the system in the vicinity of the crossover are also important because they expose fine characteristics valuable for 
understanding the system dynamics in the crossover. A powerful tool to study this is the time evolution of ground 
state survival probability [lj, [13]. It describes the dynamical behavior of the system's ground state under small 
perturbation in parameters in the system Hamiltonian. In our problem, we parameterize the Hamiltonian using the 
interaction strength modulation depth 70 such that H = if (70). If 7o has a small variation S-f so that 70 —7 70 + £7, 
the ground state's survival probablity is defined as 

M(t) = \{G( 7o )\exp{-iH^ + 5 1 )t}\G( 7o ))\ 2 . (17) 

The survival probability can be considered as a special case of the quantum Loschmidt echo 15]. Incidentally, the 
numerical method we used to simulate our system, the TEBD algorithm, is very convenient in calculating the system's 
time evolution and hence the survival probability. 

We calculated the ground state survival probability M(t) and plot the results in Fig. [7] for N = 6 and the per- 
turbation in 70 is £7 = 0.1. M(t) exhibits roughly sinusoidal oscillations in time. The amplitude of the modulation 
reaches the maximum near critical points which are consistent with those found in previous static study and shown 
in Fig. [6] Indeed, the curve for the oscillation amplitude of the ground state's survival probability can be used to 
predict the overlap between the perturbed and original ground state wave functions which is plotted in Fig.|6]Jb). The 
larger the oscillation amplitude in M(t), the more sensitive the ground state wavefunction to the perturbation in the 
Hamiltonian. 



V. SUMMARY 



In conclusion, we have made a systematic investigation of a condensate on a nonlinear ring lattice. Our studies show 
that the properties of the system are sensitive to the modulation period d of the interaction strength. In particular, the 
meanfield symmetry breaking phase transition is second order when the modulation period d is 2 but first order when 
d > 2, due to competing mechanisms present in the system driving these transitions. Our full quantum mechanical 
treatment based on the TEBD method reveals the behavior of many important quantities that are essential to the 
characterization of the system physics. 
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FIG. 7: (Color online) The amplitude of the oscillation of the survival probability M(t) for (a) d = 2 and (b) d = 3. The inset 
shows the dyanmics of M(t) for different values of 70. 

That the mean-field symmetry breaking phase transition changes from second to first order when d changes from 2 
to larger than 2 is somewhat surprising. The change of the order of the phase transition may be related to the change 
of the length scale associated with the modulation of the scattering length: For d > 2, this length scale is smaller 
compared with d — 2. Recently, Mayteevarunyoo et al. studied the symmetry breaking transition in a BEC subject 
to a nonlinear double-well potential and found that the width of the nonlinear potential plays an important role in 
controlling the transition [16j, a phenomenon that may be related to what we discovered in the current work. This is 
certainly one of the peculiar properties of nonlinear potentials that deserves further investigation. 
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Appendix A: Proof of the equivalence between the algebraic type and time-independent Gross-Pitaevskii 

equation 

Following the argument in Sec. I, we obtain the Hamiltonian in the limit N 1 as 

ii V^ 2 t t , Ho t J fx 

n = ) J —a? k a\a.kai + ^7 ) / a k a{a m an (pd+m+n-k-l 

kl klmn 

S-d+m+n-k-l) j (Al) 
which can be written in a simplified form: 

H = — aijkia\a\akai . (A2) 

ijkl 

In Sec. IIII1 we obtain the algebraic type GP equations (fl"0|) and (p~T|) which we rewrite here as: 

co = ^2 o-iju uoiUojU^u^ , (A3) 

ijkl 

= ^2 {a-ijki + ctijik) [u 0l u jUQ k (ai - 11^X0)] ■ (A4) 

ijkl 

Here, there is the unitary relation: J^i u oi u oi — 1- Based on Eqs. (|A1[) and (|A3|) . we have: 

c = 2J a-ijkiUoiU 0j ul k u* Ql = A + B, (A5) 

where A = Y,ki k2u ^u iul k u* m = J2k k2u oku'o k and 5 = 22° £] fe;mn uo k u ol ul m ul n (6 d + m + n -k-i - S-d+m+n-k-i)- 
Furthermore, by considering Eq. (|A4|) . the following relation can be derived: 



^ k 2 u 0k u Q i(uQ k ai + Uo/a fc ) + — ^ ^ u ok u a i{u* 0m a n + u* 0n a m ) (S d + m +n-k-i - &-d+m+n-k-l) = 2c xo ■ (A6) 



«7o 

klmn 



By taking advantage of the unitary relation: xo — uoiai, Eq. (|A6|) can be decomposed into the following algebraic 
equations depending on the boson operator a;: 

1 2 Uqi + ^2 u 0kU0nUQ m (5d+m+l-k-n — & -d+rn+l-k-n) = ^UOl j (A7) 

fcmro 

where // = 2c - A = c + f^iV. 

The time-independent GP equation describing Bose gases on a ring with periodic scattering length is: 

- j^W ( ) ~ 2lT ^ sin W 1^ (^)l 2 ^ (*) = ^ (*) • ( A8 ) 

with the boundary condition: -0 (#) = V' (2tt + &). By replacing (0) = -7= UQie~ lW into Eq. (|A8I) . the algebraic 
equations same as Eq. (|A7|) can be derived. This demonstrates the equivalence between Eqs. (fTCTj) and (TTTT) and the 
usual GP equation. 

Appendix B: Finding Excitation Spectrum 

The Bogoliubov spectrum of quasiparticle excitations can be determined by diagonalizing operator K as defined 
in Eq. ©. However, here, we can not obtain the representation of Eq. © using the modes jxpXmj because the 
unitary transformation U is unknown except for its matrix elements in the first row. This difficulty can be overcome 
by representing the last two terms at the r.h.s. of Eq. ((SJ) using the original mode operators |aj\a m |: 

), (Bl) 

mn inn 

where £) m „ B mn a m a n = Y.i,j,k,i a iM \ a \ a ]) '■ a k a i ■ and 
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y^m«aL a " = ^2 ai 3 kl {( a i afc ) : a ] ai : + ( a ] a; ) : a i ak '■ + ( a l a ') : a J a *; : + ( a j a *:) 1 a ! a ' : } - l»(yi a n a " _ AT) ' 



Here (A) refers to the amplitude of the operator A projecting onto the operators XoXOj XoXo? or XoXo an d : A : refers 
to the corresponding component in the operator A orthogonal to the operators (XoXo? XoXo> XoXo)- For instance, 



a\a^ : a k a t := u 0l u oj (a k - u^Xo) (a-i - u^xo) 



To find the energy spectrum for quasiparticle excitation, a generalized Bogoliubov method ( see reference 175) can 
be used to diagonalize Eq. (IBlj) . Briefly, by introducing the row and column vectors 



(B2) 



and defining the matrix 



we may rewrite the operator K as: 



M=[ ) (B>!) 



K=(co-fx)N + ^M<;-^TrA. (B4) 

Here, the matrix A = [A mn ] and the matrix B — [B mn \. Now K can be diagonalized by introducing the appropriate 
canonical transformation T: j3 = Tq. Finally, the operator K has the diagonal form: 

K = (co - n)N + ^rjTTiMT- 1 ^ - ^TrA , (B5) 

where the matrix r] = J© (— /) and Tr\MT~ x is a diagonal matrix. 
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